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O \ Abstract 

o 

■ After recent results from solar neutrino experiments and KamLAND we can def- 

initely say that neutrinos from SN1987A underwent flavor conversion, and the con- 
version effects must be taken into account in the analysis of the data. Assuming the 
normal mass hierarchy of neutrinos we calculate the permutation factors p for the 
Kamiokande-2, 1MB and Baksan detectors. The conversion inside the star leads to 
p = 0.28 — 0.32; the oscillations in the matter of the Earth give partial (and different 
CN ' for different detectors) regeneration of the original Uf. signal, reducing this factor down 

to 0.15 - 0.20 (at = 40 MeV). We study in details the influence of conversion on the 
, observed signal depending on the parameters of the original neutrino spectra. For a 

given set of these parameters, the conversion could lead to an increase of the average 
energy of the observed events up to 50% and of the number of events by a factor of 2 
, at Kamiokande-2 and by a factor of 3 - 5 at 1MB. Inversely, we find that neglecting the 

\ conversion effects can lead up to 50% error in the determination of the average energy 

' of the original Ve spectrum and about 50% error in the original luminosity. Compar- 

ing our calculations with experimental data we conclude that the Kamiokande-2 data 
Qh! alone do not favor strong conversion effect, which testifies for small difference of the 

original and spectra. In contrast, the combined analysis of the Kamiokande and 
1MB results slightly favors strong conversion effects (that is, large difference of the 
original spectra). In comparison with the no oscillation case, the latter requires lower 
^ ' average energy and higher luminosity of the original Ve flux. 
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1 Introduction 



After recent results from solar neutrino experiments, and first of all SNO PIEIISI, as well 
as from the reactor experiment KamLAND [Ij, which have selected the Large Mixing Angle 
(LMA) MSW solution of the solar neutrino problem, we can definitely say that neutrinos 
from SN1987A got converted. The neutrino flavor transformations influenced the signals 
observed in 1987 by Kamiokande-2 (K2)[Slin], 1MB [HIH] and Baksan p. Conversion effects 
must be taken into account in the analysis of the data and in the determination of the 
properties of the original neutrino fluxes. Results obtained without conversion are to some 
extent incorrect. 

The conversion of neutrinos associated to SN1987A has been extensively studied before 
(see fn]-[221 for an incomplete list of references). Here we comment on some papers which 
are relevant for our present discussion. 

It was suggested [22] that the difference of signals detected by the K2 and 1MB detectors 
can be explained partially by the oscillations of neutrinos in the matter of the Earth since 
the distances crossed by neutrinos on the way to these two detectors were different. The 
suggestion implied, however, a large lepton mixing, which was not a favored idea at that 
time. The detailed calculations have been done 13 years later jTHj, when certain hints 
appeared that LMA could be the correct solution of the solar neutrino problem. 

In connection to SN1987A, Wolfenstein considered antineutrino conversion in the star 
in the case of large mixing ^U]. He concluded that conversion leads to a harder energy 
spectrum of the observed events and, possibly, to a larger number of events. 

In the attempt to restrict the large lepton mixing, the conversion of antineutrinos in the 
non-resonance channel has been considered in details ^H]- From the analysis of the SN1987A 
data a bound on the permutation factor {p < 0.35), and consequently,on the mixing angle 
has been obtained. It was found that the bound is weaker in the LMA range, and the Earth 
matter effect further relaxes it. 

A detailed analysis of the SN1987A data based on the Poisson statistics has been per- 
formed by Jegerlehner et al. who found that the data do not allow a definite conclu- 
sion on the oscillations hypothesis. In the event that a large neutrino mixing is confirmed 
(as it has been recently), the data analysis would point toward average neutrino energies 
(at production in the star) lower than what theoretically predicted. By combining solar 
and SN1987A data, the authors of refs. [201 1^ concluded that the LMA region was the 
most suitable, among the large mixing solutions of the solar neutrino problem, to reconcile 
SN1987A data and predictions from numerical supernova codes, in agreement with [H)] . 

In ref. [IH], following the early suggestion [23], we have considered the possibility that 
certain features of the energy spectra of the events detected by K2 and 1MB can be explained 
by neutrino oscillations in the matter of the Earth. This fixes several bands in the Am^ — 
cos 26 plane. It was concluded that the data favor the parameters of LMA solution and the 
normal mass hierarchy. The inverted mass hierarchy is disfavored, unless the 1 — 3 mixing 
angle is very small fH] (see however |22j). 

The combined analysis of results from solar neutrino experiments and KamLAND lead 
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to the values of oscillation parameters 

Am^ = 7.ltli ■ 10-^eV^ tan^ 9 = 0.40 ± 0.10 , (1) 

which coincides with the third band (from the bottom in Am^ scale) found in ^H] • 

In this paper we revisit the conversion of neutrinos from SN1987A using the latest infor- 
mation on neutrino mass spectrum and mixing. We address the questions of how neutrinos 
were converted, how conversion modified the observed signals and what could be the error 
in the determination of the original neutrino fluxes if conversion is not taken into account. 

The analysis of the SN1987A data and a determination of the original spectra as precise 
as possible are needed not only to understand what happened in 1987 but also to compare 
with the results of future detections of neutrino bursts from supernovae. Detections of su- 
pernova neutrinos are rare events and furthermore each supernova is unique. Indeed, the 
mass of the progenitor, luminosity, rotation, magnetic fields, chemical composition can be 
substantially different, and, as a consequence, the properties of the neutrino fluxes can vary. 
Thus, future high statistics detections are not expected to reproduce the same features as 
those of SN1987A, but will give somehow complementary information. The comparison of 
neutrino signals from different supernovae would be extremely important for understanding 
the latest stage of star evolution, the dynamics of core collapse and explosion. 

The paper is organized as follows. In sec. 2 we consider conversion of antineutrinos and 
calculate the Ue survival probabilities and permutation factors for Kamiokande-2, 1MB and 
Baksan. In sec. 3 the effects of conversion on the observed signals are studied depending 
on the parameters of the original spectra. In sec. 4 we compare the predictions with the 
experimental results and make some indicative conclusions on the properties of the original 
fluxes. The results are summarized in sec. 5. 



2 Neutrino conversion. Permutation factor 

As a consequence of the equality of the original z/^ and Ur fluxes, = F^ = the 
conversion effects in the antineutrino channel are described by the unique probability Pee = 
P{ue i^e) which is the total Ue survival probability from the production point to the 
detector [21] (see however j2Sl)- Pee takes into account the conversion/oscillation effects 
inside the star, on the way from the star to the Earth and oscillations in the matter of the 
Earth. Using Pee? the electron antineutrino flux at the detector, Fg, can be written in terms 
of the original fluxes as 

Pg = P° + (1 - Pee)AP°, (2) 

where 

AP° = P° - P° (3) 
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is the difference of original fluxes. The combination 

P=l-Pee (4) 

is often called the permutation factor (indeed, if Pee = 0, p = 1 and Fg = that is, as a 
result of conversion the initial spectra are permuted). 

After the recent determination of the oscillation parameters, and with a few plausible 
assumptions, the physics of conversion is basically determined and the probability Pee = 
Pee(-E') can be calculated reliably. 

In general, taking into account the loss of coherence between mass eigenstates on the 
way from the star to the Earth, we can write 

^^e = 5^Per Xi^r*^ (5) 

i 

where Pgf^ is the probability of the Ug —>■ i^i transition inside the star, and p.^*^^*^ is the 
probability of z/j Vf. transition inside the Earth. Here Ui are the neutrino mass eigenstates. 
We assume the following. 

• There are only three neutrinos, and if additional (sterile) neutrinos exist, their mixings 
with active neutrinos are negligible. 

• Neutrinos have normal mass hierarchy, or if the hierarchy is inverted, the 1-3 mixing 
is very small, so that its effect can be neglected (see [21|)- 

• The density profile in the region of the 1-2 resonance coincides with the progenitor 
profile during the neutrino burst (~ 10 s). It is not affected by shock wave propagation. 
Even if the profile changes with time, the adiabaticity character of the conversion 
remains unchanged due to the large 1-2 mixing. 

(We will comment on effects of possible relaxation of these conditions in the sec. 5.) 
Under the above assumptions the dynamics of neutrino conversion is completely deter- 
mined: 

1). Inside the star the electron antineutrinos are adiabatically converted into the mass 
eigenstate vi j21j 

Ve-^Vi, (6) 

and consequently, = 1, whereas P^^ = for i 7^ 1. Indeed, in the assumption of 

normal mass hierarchy there are no level crossings in the antineutrino channel. Furthermore, 
in the production point the mixing is strongly suppressed, so that there we have Ue = i^im, 
where uim is the eigenstate of the neutrino propagation in matter. For the parameters 
the propagation is highly adiabatic and so the neutrino state coincides with during 
the whole propagation. At the surface of the star (zero matter density) equals Ui. 
Corrections due to possible deviations from the adiabaticity are negligible. 
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2) . Being an eigenstate of the Hamiltonian in vacuum, ui propagates without any change 
from the surface of the star to the Earth. 

3) . Inside the matter of the Earth z/i oscillates, in particular to Ue'- 

Vl^ Ve- (7) 

We denote the probability of this transition as Pfj^^^^. 



According to this physical picture and eq. the survival probability can be written 
immediately: 

Pee = Plf'^. (8) 

Thus, the total survival z/g '^e probability simply coincides with p^^"^*^ inside the Earth. 
If the star is not shielded by the Earth, we find 

Pi?^*"=|t/el|'=COs2^13COs2^12 (9) 

in the standard parameterization of the vacuum mixing matrix (see e.g. [SEl)- Furthermore, 
if the 1-3 mixing is small enough the probability equals 

Pee ^ cos' ^12 , (10) 

where the angle 6*12 is immediately obtained from the analysis of the solar neutrino data 
and KamLAND results in the two neutrino framework, eq. (P). 

If the 1-3 mixing is near the present upper bound, sin' 6'i3 ^ 0.06 j2Zl 1211 123 ED] , and 
we are interested in ~ 10% corrections, the complete expression Q should be used. The 
solar neutrino data allow us to measure immediately the combination cos"^ ^13 sin' 612 in high 
energy experiments and cos'^^i3(l — 0.5 sin' 2^12) in future low energy experiments. These 
combinations differ from that in eq. Q, therefore independent measurements of ^13 are 
needed to know p^j^^^^ accurately. In the present discussion of the SN1987A signal we can 
neglect corrections due to non-zero 1-3 mixing, since the errors produced in doing so are 
within the present accuracy of determination of 612. 



The probability of oscillations inside the Earth can be written as 

Pi^-*'^^ |t/,i|' + /,,„ (11) 

where freg is called the z/g regeneration factor. Notice that freg > 0, because inside the 
Earth the matter density is higher than in the cosmic space. So, there is a partial return to 
the initial condition of high density in the production point, where the neutrino state is z/g. 
Thus, the oscillations inside the Earth weaken the net conversion effect. 

In the two neutrino context, the permutation factor can be written then as 

P = sin' 6*12 - freg, (12) 
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and the flux of tlie electron antineutrinos at the detector equals: 




+ pAF' 



^F^+{smH,2-freg)AF'. 



(13) 



In spite of its simplicity, this expression gives a precise description of the neutrino conver- 
sion effect. The permutation factor does not depend on time during the neutrino burst, in 
the assumption that the matter density profile in the star does not undergo a drastic time 
evolution. 

We have performed exact numerical calculations of the probability p^^"^*^ and the regen- 
eration factor using the Earth density profile from [SJ . The fig. ^ shows the permutation 
factor p = X — p^j^^'ih Kamiokande-2, 1MB and Baksan detectors for two different choices 
of the oscillation parameters from the LMA allowed region. To remove the unobservable 
fast oscillations we have averaged the probability over energy as follows: 



In the figures Ai? = 2 MeV was taken for illustration. 
The following remarks are in order. 

Due to the larger distance crossed by neutrinos inside the Earth, for the 1MB detector 
the frequency of the oscillatory curve in the energy scale is twice as large as the frequency 
for K2. The depth of oscillations is nearly the same for all detectors, but the phase of 
oscillations is different. 

Also the averaged permutation factor p has the same energy dependence for all detectors: 
For the best fit values of oscillation parameters (upper panel) it decreases from 0.28 down 
to 0.25 at E = 25 MeV and 0.23 at E = 40 MeV. 

For K2, at low energies {E < 25 MeV), the oscillations in the matter of the Earth are 
essentially averaged out. For higher energies the averaging is not complete. Notice that the 
minima of the oscillatory curve are at E = 30 [p = 0.21 ) and 40 MeV {p = 0.19). For 1MB 
the minima are at 36 MeV and 43 MeV (p = 0.19). 

For comparison we show also the permutation factors for Am^ = 5.2 ■ 10~^ eV^ which is 
at the border of the allowed region and in the lower band identified in ^Hj from the spectral 
features. 

Notice that with the decrease of Am^ both the depth of modulation and its period 
increase. In particular, f^eg oc 1/Am^. Now in the minimum of the K2 curve at 40 MeV we 
find p = 0.15. The minima are also slightly shifted in energy. 

The permutation factor increases with the mixing according to eq. (fT2|l : at the same 
time the regeneration factor changes weakly. 

The fig. El shows values of the regeneration factor freg for K2 and 1MB in the Am^ — 
sin^ 6*12 plane for E = 40 MeV. For other energies the result can be obtained from fig. |2]by 
rescaling of Am^: freg{E,Am'^) = freg{Am'^ x (40MeV/-E)). For the best fit values of the 




(14) 
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oscillation parameters we have freg — 0.08 at Kamiokande-2 for E = 40 MeV. 
Concluding this part we can definitely say, that 

• the permutation factor due to conversion inside the star is about 0.2 - 0.4 with most 
plausible value 0.28 - 0.30; 

• for all detectors this factor is modulated significantly by the Earth matter effect, which 
suppresses the permutation down to ~ 0.17 for E = 40 MeV. The modulation effect 
due to oscillations inside the Earth increases with energy. 

The effect inside the star is determined mainly by the mixing sin^ 6*12 and practically 
does not depend on Am^. In contrast, the oscillations inside the Earth are very sensitive 
to Am^ and also depend on mixing. Future operation of the KamLAND experiment will 
allow us to determine Am^ even better, and therefore to have more precise determination of 
the Earth matter effect. Improvements in the determination of the mixing angle may follow 
from further operation of SNO. 

Notice that in spite of the relatively small values of the permutation factor and even 
smaller Earth matter effect, the influence of conversion on the neutrino fluxes can be strong. 
Indeed, if the original Ux flux has substantially larger average energy than the Ue flux, then 
Fg -C in the high energy tail of the spectrum, due to the exponential decrease of the fluxes 
with the increase of energy. Therefore the observed flux at high energy will be dominated 
by the converted flux (second term in eq. 0). 

3 Conversion effects and original spectra 

According to eq. (0) the effect of conversion on the neutrino spectra, and, consequently, 
on the spectra of observed events, is proportional to the difference of the original z/g and 
Ux fluxes. In view of the large uncertainties in these original fluxes (see e.g. the summary 
in jnSj) we will study systematically the conversion effects depending on the values of the 
parameters of the original spectra. 

We will use the following parametrization of the instantaneous original fluxes of neutri- 
nos, which depends explicitly on the integral characteristics |32j : 

where Lz is the total (integrated over the neutrino energy) luminosity in the flavor z and 
Eqz is the average energy (as it can be checked by explicit calculation); D is the distance to 
the star; and plays the role of a pinching parameter. It is estimated to be in the interval 

a^~2-5. (16) 



6 



The width of the spectrum (fT3j) equals 



The Fermi-Dirac spectrum corresponds to ^ 2.3. 

The parameters of neutrino fluxes change with time during the neutrino burst: £"02 = 
Eoz{t), Lz = Lz(t), az = az(t). In particular, the spectra may have two different time 
components (see e.g. the review in jHS] and references therein): one from the accretion 
phase, and another one from the cooling phase. During these phases the parameters of the 
spectra can be substantially different whereas within each phase they change slowly. 

The integral characteristics of the original spectra relevant for our discussion are the 
average energy of the electron antineutrino spectrum, Eoe] the luminosity in the electron 
antineutrinos, L^; the ratio of the average energies of the muon/tau and electron antineutrino 
spectra, r^, and the ratio of the corresponding luminosities, r^: 

rE^^, r,^^; (18) 
the pinchings of the Ue and Ux spectra, ae and Ox- 

In the following we consider the signals observed in the Cerenkov detectors K2 and 
1MB. The features of the predicted neutrino signal at Baksan and the consistency with the 
observed data are discussed at the end of sec. El The differential spectra of the detected 
positrons from +p — > + n are given by 



de 



N^J^"^ de'n\e, e')8\e') J dEF^iE)^^^^ , (19) 



where i = K2, 1MB, e and e' are the observed and the true energies of the positron re- 
spectively, is the number of target particles in the fiducial volume, is the detection 
efficiency. In p9|) the energy resolution function, 7?.*(e, e'), is taken in the Gaussian form 
with a, = 0.87v/e/MeV for K2 P and a, = 1.16^e/MeV for 1MB [Zj. The fluxes of 
antineutrinos at the detectors, F^{E) are given by eqs. (|TH|l and (|T!^ . We use the differen- 
tial cross section of the detection reaction, da{e' , E)/de' , calculated in [HI]. The conversion 
effects (p*) have been found using the best flt values of oscillation parameters given in eq. (^. 

The signals at the various detectors will be characterized by the following integral char- 
acteristics. 

- rate, or total number of the observed events N^: 

r^°° dN' 

N'^ (20) 
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the average energy of the detected (electron antineutrino) events, 



de e 



dN' 

~dr'' 



(21) 



- width of the observed spectrum: 



where 



(6^ 
(6-)2 



1/2 



(22) 



(e = — 



de e^ 



de 



(23) 



We calculate these characteristics of the K2 and 1MB signals as functions of the pa- 
rameters of the original neutrino spectra. The results can be considered as instantaneous 
characteristics of observed events or as the integral characteristics if the parameters of the 
original spectra do not change during the neutrino burst. This matches the constant tem- 
perature model presented in the upper panel of Table 4 in ref. ^■ 

Results for the observed numbers of events and average energies are presented in the figs. 
OHni and summarized in the Tabled They refer to the case in which ae = o-x = 3, so that 
the corresponding widths are equal and fixed to the value 7 = 0.5, according to eq. (fT7j) . 
Fig. [7| shows results for the widths F* and illustrates their dependence on a = = a^- A 
further illustration of the a-dependence is given in Table |21 



Let us summarize the effects of flavor conversion: 



1). The conversion leads to an increase of the average energy of the observed events 
(fig. El) . According to fig. El for K2 the dependence of ? on Eq^ is linear in a wide interval of 
energies (-Eoe > 5 MeV) due to the low energy threshold, and it can be approximated with 
good accuracy by 

e^ A{rE,rL)EQ, + B{rE,rL). (24) 

We find A(ri5,l) = 1.24, 1.32, 1.47, 1.58 for r^^ = 1.0, 1.2, 1.6, 2.0 respectively. For 1MB, 
due to the higher energy threshold, the the dependence is linear in a restricted interval only. 

Due to their dependence on the difference of the original Ve and Vx fluxes, eq. the 
conversion effects increase with te- In the wide energy range £^oe = (8 ~ H) MeV (and for 
tl = 1), both and e^^^^ increase by up to ~ 37% with the variation of te between 1 
and 2. This change is substantially larger than the la experimental interval, which equals 

"'^ Strictly speaking, the results in \IVA\ are not directly comparable to ours, since they were obtained using 
a different parametrization of the neutrino original fluxes. 
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Table 1: Effect of tlie neutrino conversion on characteristics of observable K2 and 1MB 
signals and extracted parameters of the original neutrino spectra. The numbers of events 
and luminosities correspond to Eq^ = 11 MeV and = 1- 

~ 7.2 — 7.5%. As shown in fig. El for Eq^ = 11 MeV the values of the average energies are 
in the intervals e^^ = 17,5 _ 23.5 and e^^^^ = 27.5 - 37.7. 

The dependence of e* on r^;, is mild at low-intermediate energies and becomes slightly 
stronger at high Eq^ and for larger te- For te = 1-6 we have that as tl varies from 0.67 
to 1.5 the change in ? does not exceed ~ 7% for K2 and ~ 3.5% for 1MB. The latter 
dependence is weaker due to the higher 1MB energy threshold, which results in a reduced 
sensitivity of the 1MB signal to the softer component of the spectrum due to the original 
flux. 

The plot (fig. El) allows us to estimate the errors in the determination of Eq^ in the case 
when conversion is not taken into account. One can see that in the no-oscillation case the 
observed average energy at K2 corresponds to i?oe = 8-8 MeV, while the 1MB observation 
is reproduced for Eq^, = 14.7 MeV. With conversion these two values become as low as 
£"06 = 5.4 MeV and Eqs = 7.6 MeV respectively. 

So, not taking into account conversion can lead to an error in the determination of the 
average energy (or temperature) of the original spectrum as large as 40 - 50%. 

2). The conversion leads in general to an increase of the number of events (an exception is 
the case of low non-electron flux, < 1, and te ~ 1). The number of events is proportional 
to the Ue luminosity, Lg-. 

N' = Nl,{Eoe,rL,rE)^, (25) 

where A'^^ is the number of events in the detector i for a fixed luminosity Lq. For definiteness 
we take Lq = 5.3 ■ 10^^ ergs, which is the integral luminosity in the Ue species determined in 
ref. inni from the analysis of data in the constant temperature scenario without oscillations. 
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In fig. m we show the dependences of Nq for K2 and 1MB on Eq^ for different values of 
the ratios and te- Also shown are the lines which correspond to the numbers of events 
observed in these detectors together with the la bands. 

While iV(f^ increases nearly linearly with Eoe (for Eq^ > 6 MeV), the dependence of 
A^^^^ is stron ger as an effect of the higher energy threshold and slow rise of the detection 
efficiency with energy. 

For Eq^ = 11 MeV and = 1 the increase of Nq with does not exceed ~ 30% for 
K2, while it can reach a factor of 3.5 at 1MB (see Table El). The increase is more sizable for 
lower Eoe, as can be seen in fig. HI 

The numbers of events increase with ri: we see a change by ~ 30% at K2 and by ~ 75% 
at 1MB for Eoe = 11 MeV and te = 1.6. 

In fig. Elwe show the dependences of the total numbers of events, A^*, on L^. for Eq^ = 11 
MeV and different values of the ratios and te- This plot allows us to evaluate the error in 
the determination of the original luminosity if the oscillation effect is not taken into account. 
The detected number of events N^"^ = 12 at K2 can be reproduced without oscillations for 
Lg = Lq = 5.3 • 10^^ ergs. With the increase of the conversion effect (increase of te or/and 
ri) the required luminosity decreases (see Table 1). E.g. we find: L^/ Lq = 0.78 for ve = 1.6. 
For 1MB the effect is much stronger: The observed N^^^ = 8 events can be obtained for 
Lg/Lo = 2.75 with no oscillations. The required luminosity decreases down to Le/Lq = 1.2 
for Te = 1.6. 

3) . The conversion effects for K2 and 1MB differ due to 

- different Earth matter effect; and 

- different experimental characteristics: thresholds and efficiencies of detection. 

In fig. El we show the ratio of the numbers of events at K2 and 1MB, = N^"^ /N^^^^ , 
as a function of Eoe. The experimental value is = 1.5. The ratio increases with the 
decrease of Eq^. The conversion leads to a decrease of r^, thus allowing to reconcile the 
required energy £"06 and the correct value of N^"^ /N^^^ . Taking Eq^. = 11 MeV and te = 1 
we find that varies between 4.2 and 1.5 depending on te (Table HI). 

In contrast to the numbers of events, the energies of the observed events at K2 and 1MB 
are modified by the conversion in rather similar ways, so that their ratio is not affected 
significantly: g-f^2^g/A/B _ g g2 _ 0.66 (see Table and fig. |31). The change is even weaker 
for lower Eq^- 

The difference of conversion effects can be used to improve the agreement of the 1MB 
and K2 signals (secEJ. 

4) . In fig. [71 we show the dependence of the widths of the observed positron spectra, F*, 
on rE for different values of the pinching parameter a = ae = olx and = 1 (results for 
other values of r^, are almost identical). It appears that this dependence is approximately 
linear for te ^ 1.2. Conversion (i.e. te > 1) leads to widening of the K2 and 1MB observed 
spectra, and the widening is larger for larger r^, reaching ~ 25 — 30% for r^; = 2. We see 
a ~ 10% decrease of T^"^ and r-'^^^'^ as ae increases from 2 to 5 (i.e., as the original spectra 
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Table 2: Effects of different values oi a = ae = dx on the observables in the no-oscillation 
case with parameters from Loredo & Lamb (lines 1-7) and in the case of oscillations 
with parameters from our concordance scenario given in sees. El and |3] (lines 8 - 14). 



become more pinched). 

The dependence of the observables on a is summarized in the Table El We show results 
for two scenarios: (1) the no-oscillation best fit scenario of ref. jB^, with Eq^ = H MeV, 
Le = 5.3-10^^ ergs and (2) the "concordance" scenario (see sec. |3)) with conversion, i?oe = 8 
MeV, Le = 8 ■ 10^^ ergs, = 1.6 and = 1. 

With the variation of a, all the quantities in the table change by at most 10%, with the 
exception of the number of 1MB events, which can vary even by a factor of 2 with respect 
to the case a = 3. In particular, the decrease of a down to a = 2, leads to an increase of 
j\fiAiB _ '^g That is, in the "concordance" scenario a smaller a, a ~ 2, could fit 

the numbers of events better. This however is compensated by the worsening of the fit to 
the average energies, so that the global fit does not improve. 

4 SN1987A data and neutrino fluxes 

Using the results of the previous sections we now compare the data with predictions and 
make qualitative statements on the interpretation of experimental results and the determi- 
nation of the original neutrino fluxes. 

Let us consider the integral characteristics of the signals. We assume that all events are 
due to Ue interactions and background. After background subtraction (according to [HIIZDj 
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we take the total numbers of events detected by K2 and 1MB as 



N^^ = 12.0 ± 3.5, N^^^^ = 8.0 ± 2.8. (26) 
The experimental value of the average energy of the detected events equals 



where the summation runs over all observed events in a given detector. We define the error 
of the average energy as 

= (28) 

where Aef is the error in the energy determination of the individual event. 
Similarly, we calculate the width: 



= Wiv^^iH^ - 1 , (29) 



and the associated error: 



Ar 



k 



\ 



(30) 



From the data we get the average energies: 

^K2 ^ j ^ j^gY -1MB ^ 3^_g ^ 2.3 MeV, (31) 

and the widths: 

r^2 ^ 0.53 ±0.1, r^^^^ = 0.23 ± 0.06 . (32) 
We show these central values and the la bands in the figs. El- El 

To compare those results with observations ()26^ l!?T|l one needs to perform an integration 
over time, taking into account the time dependence of the parameters of the original spectra 
(alternatively, one could analyze the data in short periods of time, however the statistics is 
too small). Indeed, it is expected that during the neutrino burst there is a significant change 
in the parameters of spectra with time (see e.g. [22] and references therein). However, due to 
the nearly linear dependences of the observables on the parameters, the relative conversion 
effect changes with time very weakly. 

Let us perform time averaging of the relation (j2H). If the ratios te and vl depend on 
time weakly - which is expected during the cooling phase - then the A and B coefficients 
are nearly constant and for the time averaged energy we can write 

(e) ^ A{fE,fL){Eoe) + B{fE,fL), (33) 
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where {Eq^) = J dtEQg/At. This means that the relation ()24j). and the dependences in fig. lU 
can be considered as the relation and dependences between the time averaged characteristics. 

Let us now consider numbers of events. Assuming linear dependence of A^g on Eoe, 
Nq ^ C^{rE,rL)Eoe + G'^{rE,rL) (which may not be a good approximation for 1MB), and 
weak change or and rL with time, we can write: 



where 



[C\fE,fL){Eoe)L + G\fE.fL)] ^ I dtL{t), (34) 



From this it follows that we can use the relations and figures constructed for instan- 
taneous spectra or - equivalently - the time-independent scenario, also for the case of 
time-dependence, keeping in mind that the parameters i?oe, Le, te, tl should be considered 
as some effective parameters obtained by appropriate averaging over the time of the burst. 
This still can be used to study the compatibility of signals in different detectors, but the 
real meaning of the parameters extracted is ambiguous and it depends on the detailed time 
dependence of instantaneous characteristics. In what follows we will use the same notations 
as in sec. 3, though one should keep in mind that these are effective parameters. 

1) . Analysis of the Kamiokande-2 data only. In the absence of oscillations, from fig. El 
we extract the value of the (effective) z/g average energy: 

Eoe = 8.7 ± 0.9 MeV, (la). (36) 

According to fig. |3J for this energy and = Lq the number of predicted events equals 
A'e = 7.8 ±1.5 which is about la below of observed number, eq. ()2fi|l . The exact number of 
observed events can be reproduced if Le = 1.54Lo = 8.2 ■ 10^^ ergs. 

The conversion leads to a decrease of i?oe- For instance, for te = 1.6 and r^, = 1 we find 
from fig. 121 

Eoe = 6.8 ± 0.8 MeV, (la). (37) 

According to fig. |3J this gives N^"^ = 6.3 ± 1.5 for Le = Lq. So, correspondingly, a fit 
to the observed number of events implies even higher integral luminosity: Lg = 1.9Lo = 
(10.1 ±2.4) • 10^^ ergs, which seems rather high, also keeping in mind the low average energies 
of neutrinos. 

So, the K2 data alone disfavor strong conversion effect; though within 2a even a strong 
effect {te = 1.6 ) can be easily accommodated. 

2) . Kamiokande-2 and 1MB. For 1MB, in the absence of oscillations, we find, from fig. 01 

Eoe = U.T ± 1.9 MeV, (la) (38) 
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which is more than 3a above the K2 result, eq. ()36|) . According to fig. EJfor this energy 
and Lg = Lq the number of predicted events equals N^'^^^ = 10 + 6/ — 4 which is in a very 
good agreement with the experimental result. The precise number of the 1MB events can 
be reproduced for Lg = O.SLq. 

Thus, the 1MB signal implies about 2 times higher energies of events and 2 times smaller 
integral luminosity in comparison with K2. 

Conversion leads to a decrease of the extracted energy and luminosity. For = 1.6 and 
ri we get 

£;™ = 10.3 ± 1.7 MeV, (la). (39) 

which is only 2a above the K2 result for the same te- 

With conversion effects, the agreement of the 1MB and K2 results can improve. We can 
find a "concordance" set of parameters of the original fluxes which give a better description 
of all the available data: average energies, luminosities and widths. The latter can be 
affected signiflcantly by the oscillations in the matter of the Earth. 

The ratio of the average energies extracted from the K2 and 1MB data approaches 1 
with the increase of te: ^™/^^^ = 1.7, 1.51, 1.43 for = 1.0, 1.6,2.0. Notice however, 
that the difference does change signiflcantly for r^; > 1.6 being about 2a. Furthermore, 
with the increase of r^; both the extracted energies E^J^ and Eq^^ decrease, and this makes 
it difficult to reproduce the total numbers of events, especially for 1MB. Considering this, 
we find that the optimal energy is in the interval E^e = (7 — 9) MeV, where one can get 
la deviation from both K2 and 1MB extracted values. This result does not depend on the 
luminosity Lg. 

Considering the numbers of events, notice that once the correct ratio r^r = N^'^ /N^^'^^ 
is reproduced, the numbers of events can be obtained by fitting Lg. From the data we find 

^N^J^ = l-5^oi- (40) 

However, according to fig. IHl for the small Eq^. required by the fit of the average energies 
the ratio turns out to be too large. Indeed, for te = 1 {i.e., no oscillations), the 
"optimal" energy is Eq^ = 10.3 MeV; for these parameters the ratio of number of events 
equals tn = 5.3. We find Eq^ ~ 9.6 MeV and rjv = 5.2 for rE = 1.2; Eoe ~ 7.7 MeV and 
rjy = 5.2 for te = 1.6; Eq^ ~ 6.2 MeV and = 4.9 for ve = 2.0. So, r^v is substantially 
larger than the experimental result and does not change practically by conversion once the 
parameters are optimized by the average observed energies. 

Let us now discuss the widths of the observed spectra F*, shown in fig. U\ (for E^e = 11 
MeV and tl = 1). One can see (see also eq. (jS2I)), that the 1MB spectrum has smaller 
width than the one of K2. The former is better reproduced without oscillations (or small 
ve) and with large a, a ~ 5, while the latter points toward the case of oscillations with 
ve ~ 1-6 and a ~ 2. The widths F* decrease slightly (by ~ 10%) with the decrease of Eoe] 
we refer to Table El for results at i?oe = 8 MeV (concordance scenario). 
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To get an idea about the preferable values of parameters of the original spectra we 
calculate defined as 



for various sets of the parameters. Here = ]\[K2^ jsjimb ^^K2^^imb ^y^'^^T^^^^ ^ are the 



numbers of events, average energies and widths of the positron spectra observed at Kamiokande2 



Without oscillations and Eq^ = 11 MeV, L = Lq, = a^: = 3 we find = 16.2. 
With oscillations for E^e = 8 MeV, L = 8 ■ 10^^ ergs, = ax = 3, te = 1.6, vl = 1, 
we obtain = H-O (see Table El for the dependence of results on the pinching parameters 
a^)- We checked that the improvement in the is mainly due to the smaller deviation of 
the predicted e^^ from the observed value. This deviation amounts to about 2.6 a for no- 
oscillations, and to ~ 1.4 a for the case with oscillations. Correspondingly, the contribution 
of e^^ to the decreases from 7 to 2 when oscillations are included. Notice however that 
with the decrease of Eqe implied by conversion effects, substantially larger integral luminosity 
is required to reproduce the numbers of events. The combination of smaller average energy 
and larger luminosity corresponds to a larger radius of the neutrinosphere: Rns oc Eq^L^^'^. 
It follows that in our concordance scenario Rns is about 2.4 times larger than in the no- 
oscillation scenario, which gives Rns — 20 — 30 Km |HH|- This increase is within the range 
allowed by the several uncertainties of astrophysical nature and by individual variations 
from one star to another. 

Thus we conclude that conversion leads to a certain improvement of the global fit of 
the data. The improvement requires lower average energies of the original i^e spectrum and 
higher i^e integral luminosity. Conversion effects lead to a better agreement of the observed 
average energies with the data sample, but this improvement is compensated partially by a 
worsening in reproducing the numbers of events. 

3) . Baksan results. If all the 5 events recorded by Baksan are attributed to signal, anal- 
yses of these data without oscillations lead to neutrino luminosities which are much larger 
than those implied by the K2 and 1MB data samples (see the account in ref. [22] )• This 
indicates that part of the Baksan events are due to background, and, in absence of a specific 
prescription for background subtraction, it is difficult to include the Baksan results in our 
analysis. In ref. [22] a method to account for a background origin of part of the Baksan 
events is proposed. Notice that in presence of oscillations the K2 and 1MB data require 
larger neutrino luminosities, therefore it is possible that including neutrino conversion will 
reduce the number of the Baksan events attributed to background. 

4) . Spectrum of events. As we have marked in the introduction, the oscillation parame- 
ters we use are in the band determined previously from the fit of the energy spectra of K2 
and 1MB ^H] (see also fig. 12]). So the oscillation parameters, and in particular Am^, are 
already optimized to get the best description of the spectra of events. We mark that, in 
contrast with what discussed in fH], here the numbers of events, widths and average ener- 




2 



(41) 



and 1MB. 
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gies of the observed positrons depend very weakly (about ~ 2 — 5% at most) on the phase of 
oscillations in the Earth. This is motivated by the larger value of Am? used here, for which 
the permutation factor has faster oscillations with smaller depth. These oscillations can be 
seen in figure El (in which no averaging has been applied), however in the observed signal 
they are almost completely averaged out as an effect of the experimental energy resolution. 

In fig. IHl we show the calculated spectra of the events at Kamiokande-2 and 1MB for 
the "concordance" set of parameters of the original fluxes: -Eoe = 8 MeV, Le = 8 ■ 10^^ 
ergs, te = 1.6, = 1 (solid lines). Also shown are the spectra without oscillation for 
the best fit parameter set determined in [SB] ( "constant temperature model" ) (dotted) and 
the spectrum for the concordance set of spectral parameters without oscillations (dashed). 
Comparing solid and dashed lines we conclude that the conversion effect is significant in the 
whole energy range for 1MB and for e > 15 MeV at K2. 

At the same time the spectra with best fit values of parameters (without conversion) 
and with the concordance set of parameters and conversion are rather close. The difference 
is mainly in the low energy K2 spectrum. We stress that, though similar, the spectra with 
and without conversion require substantially different parameters of the original spectra. 

5 Discussion and Conclusions 

1) . After the identification of the solution of the solar neutrino problem and KamLAND 
results, we can definitely say that neutrinos from SN1987A underwent flavor conversion 
inside the star and oscillations in the matter of the Earth. 

In the assumption of the normal mass hierarchy, the conversion probabilities can be cal- 
culated with good precision. We find that the permutation factor is about p = 0.28 — 0.32 
due to conversion inside the star and oscillations in the matter of the Earth suppress the 
permutation. The Earth effect increases with energy, and at ~ 40 MeV p decreases down 
to 0.15- 0.20. 

2) . The conversion effects on observables depend strongly on the properties of the original 
neutrino fluxes, in particular, on the average energies and widths of the original spectra of Ue 
and i'^/t- For a given set of these parameters, conversion leads to an increase of the number 
and of the average energy of the observed events as well as of the widths of the observed 
spectra. The conversion effects are different for K2 and 1MB. By varying parameters in the 
ranges allowed by astrophysics, we have found that the average energy of events can increase 
by 30 - 50 % and the number of events by 50 % for K2 and by up to a factor of 3 for 1MB. 

Vice versa, conversion changes the average energies and luminosities of the original neu- 
trino fluxes extracted from the observations. In particular, it leads to a decrease of the 
original energy and (for fixed Eq^) of the luminosity of z/g. We find that the decrease can 
be up to ~ 50 — 70%. This leads to uncertainty in the determination of parameters of the 
original spectra. 
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3) . Comparing calculations with the real signals from SN1987A we find that the K2 data 
alone do not show significant conversion effects, and can be well described by original fluxes 
with small conversion effects. This would testify for small difference of the original Ve 
and v^j, spectra. At the same time the K2 data do not exclude strong conversion. In this 
case, however, the original spectra should have lower energies and higher luminosities in 
comparison with the no oscillation case. 

The characteristics of the original spectra extracted from the K2 and 1MB data exhibit 
substantial differences. These can be partially reduced by conversion effects: that is, the 
conversion improves the combined flt of the K2 and 1MB data. The improvement is not 
dramatic, though, and does not lead to substantially more coherent overall picture. It 
requires strongly different original spectra, Te = 1-5 — 2, low average energy of the z/g- 
spectrum, Eq^ = 6.5 — 8.5 MeV, and high luminosity: Lg = (8 — 12) ■ 10^^ ergs. 

The oscillations in the matter of the Earth substantially modify the high energy parts 
of the spectra of events at K2 and 1MB. 

Our conclusions are valid for normal mass hierarchy or inverted hierarchy provided that 
^13 is negligibly small. If the hierarchy is inverted and sin^ ^13 ^ few ■ 10~^, the u^^i^x 
permutation in the star is complete [2^ EH] , so that the z/g flux at Earth is entirely due to 
the original / Vj. flux. It follows that the parameters extracted from the data analysis 
refer to the non-electron flavors produced inside the star and therefore would lead to a 
completely different test of supernova theory, with respect to the non-oscillation case. 

In this same scenario of mixing and hierarchy, the amount of Ve^v^ permutation could 
change at late times, as the supernova shock-wave reaches the external layers of the star, thus 
modifying the adiabaticity character of the 6'13-induced MSW resonance j2Zl EH EHl EHl • The 
effect on the time integrated neutrino signal is however small |38| l39] and negligible with re- 
spect to the large uncertainties of statistical and astrophysical nature on the SN1987A data. 

4) . In general, one needs to perform an analysis which employs the whole information con- 
tained in the data: both integral and differential (energy spectrum, arrival time) as well as 
errors in the determination of the energies of events, background and angular information. 
In view of the small number of detected events, the optimal type of analysis would be along 
the lines of the work of Loredo and Lamb [SBj, with the conversion effects taken into account. 
The present study will allow to better understand the results of such a global flt. Our con- 
clusions concerning the interpretation of the observed signals have qualitative character only. 

5) . It would be important to compare these results on SN1987A with those of future SN 
neutrino detections. The latter will have high statistics and therefore will provide the 
possibility to disentangle the oscillation effects and the properties of the original fluxes. We 
mark that many of the results presented here have general character and therefore will apply 
to future data as well. In particular, the effects of conversion in the star will be valid. At 
the level of probabilities (permutation factor), the effects of oscillations in the Earth will 
probably be different, due to different trajectory of the neutrinos in the Earth. However the 
difference with respect to SN1987A may be negligible in the observed energy spectra if the 
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energy resolution of the detector is larger than the size of the spectral modulations due to 
oscillations. The analysis of future supernova data will result in a better understanding of 
the generation of neutrino fluxes and therefore in more reliable predictions of fluxes also for 
SN1987A. In this way a more precise interpretation of the SN1987A signals can be done. 
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Note added 

It has been proposed recently that the signal detected by the LSD detector at Mont Blanc 
|4()j . 4.7 hours before the Kamioka-2 and 1MB signals, is due to the primary collapse of a 
fastly rotating Fe-O-C stellar core that leads to the formation of a collapsar j3T]. The latter 
undergoes fragmentation into a binary of neutron stars, which loose angular momentum by 
emission of gravitational waves. This model predicts an early neutrino burst, produced in 
the URCA processes (neutronization), and a second one, after several hours, which results 
from the collapse of the heavier star of the binary system. 

The early burst is composed dominantly of u^s with average energies of 30 - 40 MeV. 
The signal in LSD is mainly due to CC interactions of these u^s with nuclei of Fe in the 
walls of the detector. The produced electrons generate electromagnetic showers which loose 
energy in the walls and appear in the scintillator as signals with energies of 7 - 10 MeV. 

Neutrino conversion can substantially change this picture. In the most plausible scenario 
of normal hierarchy and tan^ 613 ^ lO"'^ the z/g flux is almost completely converted to a u^^/ 
flux in the star [211 EE]- Furthermore, the regeneration effect in the Earth is negligible. So, 
the proposed signal due to the CC- interactions in the LSD detector should not appear. 

In the case of very small 1-3 mixing, tan^ ^13 <C 10~^ or inverted mass hierarchy, the 
signal is suppressed by the factor sin^^i2 ~ 0.3 (though some regeneration effect in the 
matter of the Earth is present). In this case the explanation of the LSD signal would 
require an original neutrino luminosity about 3 times larger than that in |41j . 

We mark that the LDS signal could be explained by NC interactions induced by the 
v^/fr flux. However, due to the small cross section, to reproduce the observed number of 
events a much larger original luminosity would be necessary. 
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Figure 1: The permutation factor, p = 1 — Pie, as a function of the neutrino energy at 
Kamiokande-2, 1MB and Baksan. The two panels refer to two choices of the oscillation 
parameters from the LMA allowed region, one of which (upper panel) is the best fit set of 
values in eq. (0). For illustration purposes, we have averaged p over energy according to eq. 
dUl), with AE = 2 MeV. 
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Figure 2: Regions of constant values of the regeneration factor freg at K2 and 1MB in 
the plane Am^/(10~^ eV^) — sin^ ^ (vertical and horizontal axes respectively), for neutrino 
energy E = AO MeV. The white squares denote the LMA best fit point, eq. (^. 
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Figure 3: The average energy of positrons in the detectors K2 and 1MB, as a function 
of the average energy of the original Vf. flux, Eq^., for the best fit LMA parameters (eq. (P), 
see also fig. and different values of = Eo^/Eoe (see legend) and = Lx/Lg. The 
solid, dashed and dotted-dashed lines correspond to = 1, 1.5,0.667 respectively. The no 
oscillation case (or, equivalently, the case te = 1) is shown for comparison. The dotted 
(horizontal) lines represent the values of e* extracted from the data for each experiment, 
with the la error (eq. ^^)- 
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Figure 4: The predicted numbers of events at K2 and 1MB as a function of the average energy 
of the original v^. flux, E'oe- The line identification, as well as the oscillation parameters, are 
the same as in fig. El The dotted (horizontal) lines represent the experimental result with 
the lo" error (eq. (jSH)))- We use Lg = Lq = 5.3 ■ 10^^ ergs. 
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Figure 5: The predicted numbers of events at K2 and 1MB as a function of the ratio 
Lg/Lo {Lq = 5.3 ■ 10^^ ergs) for Eq^ = 11 MeV. The hne identification, as well as the 
oscillation parameters, are the same as in fig. El The dotted (horizontal) lines represent the 
experimental result with the la error (eq. ()2fj|l ). 
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Figure 6: The ratio of the numbers of events at K2 and 1MB, r^r, as a a function of the 
average energy of the original flux, Eq^.. The line identification, as well as the oscil- 
lation parameters, are the same as in fig. 01 The dotted (horizontal) lines represent the 
experimental result with the la error (eq. (001) )• 
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Figure 7: The predicted widths of the K2 and 1MB energy spectra as a function of for 
different values of the pinching parameter ag. The oscillation parameters are the same as 
in fig. El We have taken Eq^. = 11 MeV, = 1 and = oie- The dotted (horizontal) lines 
represent the experimental result with the la error, eq. 
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Figure 8: The predicted positron energy spectra at K2 and 1MB with and without oscilla- 
tions. The solid lines correspond to E'oe = 8 MeV, Lg = 8 ■ 10^^ ergs, = 1.6, = 1 and 
the best-fit LMA parameters (eq. (^, see also fig. Q). The dashed lines correspond to the 
no-oscillation case with the same parameters: Eqs = 8 MeV, and Lg = 8 ■ 10^^ ergs. The 
dotted lines correspond to no-oscillations with the best fit parameters of ref. Eq^ = H 
MeV, Lg = 5.3 ■ 10^^ ergs. The histograms show the experimental results. 
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